Gallium adsorption on (0001) GaN surfaces 
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We study the adsorption behavior of Ga on (0001) GaN surfaces combining experimental specular 
reflection high-energy electron diffraction with theoretical investigations in the framework of a kinetic 
model for adsorption and ab initio calculations of energy parameters. Based on the experimental 
results we find that, for substrate temperatures and Ga fluxes typically used in molecular-beam 
epitaxy of GaN, finite equilibrium Ga surface coverages can be obtained. The measurement of a 
Ga/GaN adsorption isotherm allows the quantification of the equilibrium Ga surface coverage as a 
function of the impinging Ga flux. In particular, we show that a large range of Ga fluxes exists, 
where 2.5 ±0.2 monolayers (in terms of the GaN surface site density) of Ga are adsorbed on the GaN 
surface. We further demonstrate that the structure of this adsorbed Ga film is in good agreement 
with the laterally-contracted Ga bilayer model predicted to be most stable for strongly Ga-rich 
surfaces [J. E. Northrup et al., Phys. Rev. B 61, 9932 (2000)]. For lower Ga fluxes, a discontinuous 
transition to Ga monolayer equilibrium coverage is found, followed by a continuous decrease towards 
zero coverage; for higher Ga fluxes, Ga droplet formation is found, similar to what has been observed 
during Ga-rich GaN growth. The boundary fluxes limiting the region of 2.5 monolayers equilibrium 
Ga adsorption have been measured as a function of the GaN substrate temperature giving rise to a 
Ga/GaN adsorption phase diagram. The temperature dependence is discussed within an ab initio 
based growth model for adsorption taking into account the nucleation of Ga clusters. This model 
consistently explains recent contradictory results of the activation energy describing the critical Ga 
flux for the onset of Ga droplet formation during Ga-rich GaN growth [B. Heying et al,, J. Appl. 
Phys. 88, 1855 (2000); C. Adelmann et al, J. Appl. Phys. 91, 9638 (2002).]. 
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INTRODUCTION 



Due to its importance as a base material for the fab- 
rication of optoelectronic devices in the blue and ultra- 
violet spectral region, GaN has been extensively stud- 
ied in recent years. The comprehensive study of ma- 
terial and device properties^ has recently been com- 
plemented by an increasing number of studies con- 
cerning surface structures^^iSiLSiSiiSi and associated 
growth mechanisms, in particular for growth by plasma- 
assisted molecular-beam epitaxy (PAMBE) on the Ga- 
polar (0001) surface, n i 12 i 13 i 14 i 15 i 16 i 17 

One of the fundamental results of these studies is that 
GaN growth by PAMBE ought to be carried out under 
Ga-rich conditions in order to obtain a smooth surface 
morphology and optimized material properties. At low 
growth temperatures, however, this leads to Ga accu- 
mulation and droplet formation, which is detrimental 
to the GaN epilayer quality>i2iiii As a consequence, it 
was thought that optimum GaN growth conditions must 
be as close as possible to Ga/N stoichiometry. It has 
been recently observed that such Ga accumulation can 
be prevented when growing GaN at high temperatures 
and small Ga excess fluxes J2iAl In particular, it has been 
shown that, at high growth temperatures, a wide range 
of Ga fluxes exists, for which a finite amount of excess 
Ga is present on the GaN surface whose quantity is in- 



dependent of the value of the Ga fiux»i£ Such conditions 
may provide a "growth window" for GaN PAMBE, i.e. 
a region, where the growth mechanisms and the surface 
morphology are independent of fluctuations of Ga flux 
and growth temperature^ 4 *^ 

However, the quantitative description of a GaN 
"growth diagram", which describes the Ga surface cov- 
erage during growth as a function of Ga flux and growth 
temperature, has not yet been achieved: the results on 
the temperature dependence of the critical excess Ga flux 
at the onset of Ga droplet formation are contradictory, 
yielding activation energies of 2.8 eV (Ref. fl4f) and 4.8 eV 
fRef.ll7j). respectively. This suggests that the underlying 
mechanisms of Ga accumulation are not yet understood. 

To address this discrepancy, we have performed Ga 
adsorption measurements on (0001) GaN. We discuss the 
results in the framework of a lattice-gas growth model for 
adsorption, which is based on ab initio calculated param- 
eters. This model explains the origin of the apparently 
contradicting parameters derived from previous experi- 
mental studies 1 ^ 7 - and gives a consistent description. 



II. EXPERIMENTAL PROCEDURE 

The adsorption experiments were performed in a 
MECA2000 molecular-beam epitaxy chamber equipped 
with a standard effusion cell for Ga evaporation. The 
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chamber also contains a rf plasma cell to provide active 
nitrogen for GaN growth. The pseudo-substrates used 
were 2/j,m thick (0001) (Ga-polarity) GaN layers grown 
by MOCVD on sapphire. The substrate temperature T$ 
was measured by a thermocouple in mechanical contact 
to the backside of the molybdenum sample holder and 
shielded from direct heating. Prior to all experiments, a 
100 nm thick GaN layer was grown under Ga-rich condi- 
tions on the pseudosubstrates to remove the influence of 
a possible surface contamination layer. 

Ga fluxes $ have been calibrated to Ga effusion cell 
temperatures by reflection high-energy electron diffrac- 
tion (RHEED) intensity oscillations during N-rich GaN 
growth at a substrate temperature of Ts — 620 °C. In 
these conditions, the growth rate is actually proportional 
to the impinging Ga flux. It is sound to assume that 
the Ga adatom sticking coefficient is unity at such a low 
substrate temperature, which permits an absolute cali- 
bration of the Ga flux in terms of the GaN surface site 
density. 

The Ga surface coverage was assessed by analyzing the 
specularly-reflected RHEED intensity by a method de- 
scribed in Refs.lT^ andl2H This method uses the oscilla- 
tory transients in specular RHEED intensity, which are 
observed during Ga adsorption and desorption on/from 
(0001) GaN surfaces. It has been shown that the du- 
ration of these transients can be qualitatively related to 
the amount of adsorbing or desorbing Gai^^ In gen- 
eral, the relation between intensity and Ga coverage is 
unknown. Although tempting, the interpretation of these 
electron reflectivity transients in terms of RHEED oscil- 
lations is not obvious. Furthermore, it must been noted 
that the shape of the transients (albeit not their dura- 
tion) depends on the diffraction conditions, notably the 
incidence angle. Of course, the modeling of electron re- 
flection would allow to directly relate the RHEED inten- 
sity to the Ga coverage (and the surface structure) but 
this is beyond the scope of this work. 

However, the total duration of the transients occur- 
ring during Ga desorption can be used to qualitatively 
estimate the amount of Ga adsorbed on the surface. In 
Ref.GI an indirect method has been used to draw limited 
quantitative information from the desorption transients; 
below we will demonstrate a fully quantitative calibra- 
tion relating the Ga desorption transient time to abso- 
lute Ga surface coverages. This allows us to circumvent 
the problem that the RHEED intensity cannot in general 
be easily related to adatom coverage. The experimental 
procedure is thus as follows: to assess the Ga quantity 
present after Ga adsorption for different impinging Ga 
fluxes, the Ga flux has been interrupted after a fixed ad- 
sorption time. The subsequent variation of the specular 
RHEED intensity due to Ga evaporation under vacuum 
has been recorded. 
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FIG. 1: Specular RHEED intensity during Ga desorp- 
tion from a (0001) GaN surface. Beforehand, Ga adsorp- 
tion has been carried out at Ga fluxes <£> as indicated. For 
$ < 0.68 ML/s, the desorption transients correspond to equi- 
librium adsorption, i.e. they do not change as a function of the 
previous adsorption time. For $ > 0.68 ML/s, this does not 
hold and the desorption transients depend on the adsorption 
time (here 1 min). The substrate temperature is Ts = 740 °C. 



III. RESULTS 

A. Ga adsorption 

Figure ^ shows the variation of the specular RHEED 
intensity after lmin of Ga adsorption (substrate tem- 
perature Ts — 740° C, Ga fluxes $ as indicated) and 
subsequent interruption of the Ga flux (at t = 0). We 
observe oscillatory transients during Ga desorption from 
the (0001) GaN surface. To quantify the desorption pro- 
cess, we define a desorption time tdes as the time interval 
between the shuttering of the Ga flux and the last inflec- 
tion point in RHEED intensity (Fig. [TJ. The desorption 
time depends on the amount of Ga present on the GaN 
surface after adsorption. 16 The fundamental finding is 
that for Ga fluxes below <E> = 0.68 monolayers (ML)/s, 
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FIG. 2: Ga desorption time as a function of the amount FIG. 4: Ga adsorption phase diagram indicating the Ga sur- 
of nominally impinged Ga, 9, for two different Ga fluxes, as face cove r a ge as a function of impinging Ga flux $ and sub- 
indicated. Ts — 740 C. strate temperature Ts- The definition of regions 1-3 follows 

Fig. El The inset shows an Arrhenius plot of the data. 



the desorption transients — and thus the desorption time 
ides — are independent of the previous adsorption time, 
i.e. of the amount of nominally impinged Ga. This is 
consistent with the results in Ref . [l6| and is visualized in 
Fig. |21 which shows td cs as a function of the amount of 
nominally impinged Ga 9 (defined as the product of the 
Ga flux and the adsorption time) . The derivative of this 
curve gives the (coverage and flux dependent) Ga stick- 
ing coefficient. Data are shown in Fig. |2 for a substrate 
temperature of Ts — 740 °C and two different Ga fluxes: 
for <f> = 0.30 ML/s, the desorption time (and hence the 
amount of adsorbed Ga) monotonously increases until it 
saturates after 9 ~ 4 ML. At larger 9, the adsorption has 
reached equilibrium and the coverage remains constant. 

For a higher Ga flux of $ = 0.90 ML/s, we observe the 
same behavior for 9 < 4 ML, but thereafter, we find no 
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FIG. 3: Ga desorption time as a function of impinging Ga 
flux $ after equilibrium has been attained for regions 1 and 2, 
and after 1 min of Ga adsorption in region 3. The substrate 
temperature is Ts = 740 °C. 



saturation. Instead, a continuous increase of the desorp- 
tion time is observed, corresponding to Ga accumulation 
on the GaN surface. A more detailed analysis finds that 
the desorption transients for <E> < 0.72 ML/s correspond 
to finite equilibrium Ga surface coverages. For higher Ga 
fluxes, no finite equilibrium coverages exist and Ga will 
infinitely grow and finally form macroscopic droplets on 
the surface. 

These results of Fig-^are summarized in Fig.|3 which 
shows the Ga desorption time — related to the amount 
of adsorbed Ga — as a function of the impinging Ga flux 
at a constant substrate temperature Ts = 740° C. It can 
be regarded as a Ga adsorption isotherm. We can dis- 
criminate three different regions: (1) an S-shaped increase 
of the Ga coverage for $ < 0.20 ML/s, (2) a constant Ga 
coverage up to $ = 0.72 ML/s, independent of the Ga 
flux, and (3) Ga accumulation and no finite equilibrium 
Ga coverages for higher $. It is worth noting that the 
transitions between the three regimes are discontinuous 
within the experimental precision of 1 °C of the Ga effu- 
sion cell (~7x 10~ 3 ML/s for fluxes around 0.5 ML/s). 
In particular, no intermediate equilibrium coverages have 
been found between regimes 1 and 2. Therefore, the 
transition fluxes between the different regimes are well 
defined. 

To fully assess the thermodynamics of the adsorption 
process, the variation of the adsorption isotherm has to 
be known as a function of substrate temperature. In the 
following we will restrict ourselves to the study of the 
variation of the transition fluxes between the different 
regimes. The result is plotted in Fig. 01 We see that 
the transition fluxes vary exponentially with substrate 
temperature. The inset shows an Arrhenius plot of the 

(12) 

data, yielding activation energies of E y A = 5.2 ± 0.1 eV 

and E% 3) = 5.1 ± 0.05 eV for the 1 -> 2 and 2 -> 3 
transition, respectively. These values, the prefactors, and 
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TABLE I: Experimental activation energies E^ p and prefactors vf^ for the transition fluxes between different Ga coverage 
regimes as obtained from the adsorption (this work, Fig. [1J and growth phase diagrams (from Refs. as indicated), v denotes 
the GaN growth rate (the N flux), a the linear temperature coefficient of the adsorption energy, and i/^™ the renormalized 
prefactor (see Sec. IV). 



Reference 


Transition 


ET (eV) 


^ xp (Hz) 

dcs V / 


v (ML/s) 


a (meV/K) 


*S(Hz) 


This work 


1 -» 2 


5.2 


3 x 10 25 





-2.3 


4.7 x 10 25 


This work 


2^3 


5.1 


2 x 10 25 





-2.2 


1.5 x 10 25 


Ref. ,17 


1 -f 2 


3.7 


5 x 10 17 


0.28 


-0.8 


1.6 x 10 18 


Ref. F7. 


2^3 


4.8 


1 x 10 24 


0.28 


-1.9 


4.8 x 10 23 


Ref. 14 


2^3 


2.8 


1 x 10 14 


1.1 


> -0.1 


5.5 x 10 13 



the corresponding values for the growth phase diagram interval will lead to a constant offset ides, which does not 

contribute to c eq in Eq. J2J). 

Using Eq. J2J), we can now calibrate the Ga adsorp- 
tion isotherm in Fig. [3] Since the relation $des(c cq ) — 
— $(tdos) implies steady-state conditions, only the data 
in regions 1 and 2 can be treated. Applying Eq. 
to the data in region 3 does not lead to the amount of 
Ga adsorbed after a finite time but to incorrect values 
because under these conditions |$dcs| < \&\- Since ex- 
perimental data are available only for <E> > 0.032 ML/s, 
the isotherm has been extrapolated by an exponential 
for smaller fluxes. The contribution of the interval 
< $ < 0.032 ML/s to the overall integral in Eq. © 
is 0.04 ML. This is only a minor correction, suggesting 
that the specific form of the extrapolation function is 
not important. The result of the calibration is plotted 
in Fig. [21 We observe that in region 1, the Ga coverage 
increases from almost zero to a value of 0.98 ML, close 
to 1 ML. The coverage then increases abruptly to a value 
of 2.5 ML in region 2. Typical systematical errors can be 
estimated to be of the order of ±0.2 ML. 

What is the detailed structure of such adsorbed Ga 



responi 

in Refs. [IT] (1 — * 2 corresponds to B — ► C and 2^3 
corresponds to C — ► D) and Q (discussing the 2^3 
transition only) are summarized in Tab. [3 and will be 
discussed in Sec. IV. 



B. GaN surface structures 

The adsorption isotherm in Fig. |3] is given in terms of 
the Ga desorption time, which is only qualitatively re- 
lated to the amount of adsorbed Ga. For a fully quan- 
titative treatment of Ga adsorption, one must know the 
dependence of the Ga desorption rate <J>d cs on the Ga 
surface coverage c. This would allow the modeling of the 
Ga re-evaporation, as the desorption rate is given by 



$des(c) 



dc 
dt 



(1) 



with the initial condition c(i = — ides) = c oq , which de- 
notes the amount of Ga adsorbed in equilibrium condi- 
tions, i.e., before Ga desorption sets in. After the time 
interval ides, the Ga coverage becomes zero. The knowl- 
edge of ides would thus allow us to compute c oq if $des(c) 
was known (which it is not). This requirement can be 
circumvented by considering that in equilibrium, the im- 
pinging Ga flux $ must exactly balance the evapora- 
tion rate, hence < E > dcs(c cq ) = — < J 5 (c cq [idos])- Integrating 
Eq. l|T]l. taking the first derivative with respect to ides, 
and using the above substitution leads to 



'*) = 



<9<f>' 



(2) 



This expression allows the computation of c eq from id es 
as a function of $ (which is known from the experimental 
data) and can be evaluated numerically. Note that c eq 
only depends on the derivative of ides(^), which means 
that ides does not necessarily have to denote the very end 
of Ga adsorption but can be taken after any time interval, 
as long as it is well-defined in the RHEED signal. Any 
further Ga desorption after the end of the chosen time 
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FIG. 5: Calibrated Ga adsorption isotherm at Ts = 740 °C 
using Eq. The data are derived from Fig. [2] As the 

method can only be applied to equilibrium coverages, only 
regions 1 and 2 are represented. 
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FIG. 6: (a) RHEED pattern (azimuth (1120)) of a (0001) 
GaN surface with a Ga bilayer present at Ts = 600 °C. The 
white arrows indicate additional streaks due to the Ga film, 
(b) RHEED intensity profile in the (1010) direction evidenc- 
ing the additional RHEED streaks due to the Ga bilayer. 



The formation of a laterally-contracted bilayer struc- 
ture in regime 2 is further corroborated by the obser- 
vation of supplementary streaks in the RHEED pattern 
after Ga adsorption in regime 2 and rapid quenching of 
the sample down to substrate temperatures below about 
Ts = 650 °C [see Fig. EJa)]. At higher temperatures, the 
supplementary streaks are too weak to be detected unam- 
biguously. The lattice parameter corresponding to these 
streaks is found to be 2.73 ±0.03 A by a fit using pseudo- 
Voigt functions and assuming the bulk lattice parameter 
for the GaN layer. This value compares very favorably 
to that obtained by the ab inito calculations (2.75 
Combined with the adsorption results, this strongly sug- 
gests that adsorption in region 2 leads to the formation 
of a laterally-contracted Ga bilayer. 

The shape of the oscillation transients in Fig. ^ sug- 
gests that, in region 3, Ga droplets are formed on top 
of this Ga bilayer. This Ga de-wetting transition may 
indicate that the attractive interaction energy of a Ga 
adatom with the surface is maximum in the second layer 
and lower in the third layer iSi This conclusion is consis- 
tent with the first principles results which will be dis- 
cussed in the next section. Ga thus grows in a Stranski- 
Krastanow mode on (0001) GaN. 

Finally, the different adsorption regimes can be sum- 
marized as follows: (1) Ga coverage c < 1 ML, i.e. suc- 
cessive formation of a coherent Ga monolayer, (2) a Ga 
coverage of c = 2.5 ML, forming a laterally-contracted Ga 
bilayer, and (3) Ga accumulation and droplet formation 
on top of a Ga bilayer. In the following, we will derive an 
ab initio based growth model which describes the tem- 
perature dependence of the transition fluxes between the 
different regimes. 



IV. DISCUSSION 



films? It must be kept in mind that the Ga fluxes have 
been calibrated by GaN RHEED oscillations and are 
hence given in terms of the GaN surface site density. 
A GaN coverage of 1 ML thus indicates the adsorption 
of a Ga adatom on each GaN site. This suggests that 
Ga adsorbs in region 1 as a coherent (pseudomorphic) 
adlayer. 

At the 1 — * 2 transition, the Ga coverage increases 
by about 1.5 ML, i.e. by more than a pseudomor- 
phic adlayer. This compares favorably to the laterally- 
contracted bilayer model, which has been calculated to 
be the most stable structure of Ga-rich (0001) GaN 
surfaces^ It consists of two adsorbed Ga adlayers on top 
of a Ga-terminated GaN surface. The first layer is found 
to be pseudomorphic to the GaN surface but the second 
one has an in-plane lattice parameter of 2.75 A, about 
13.8% smaller than that of GaN (3. 189 A). The second 
layer thus contains 1.3 ML of Ga in terms of the GaN 
surface site density, in good agreement with the experi- 
mental value of 1.5 ML. 



An intuitive connection between adsorption and 
growth phase diagrams can be made by assuming that 
the excess Ga in Ga-rich growth conditions behaves as if 
it would be adsorbed on a GaN surface. Then, adding 
the GaN growth rate (0.28 ML/s in the experiments of 
Ref. ITt^ to the adsorption phase transition fluxes should 
reproduce the growth phase diagram of Ref. 17. This 
comparison is shown in Fig. \7\ It demonstrates good 
overall agreement, suggesting that the assumption is 
valid. However, it must be noted that the critical fluxes 
derived from adsorption measurements fall below those 
in the growth phase diagram, although differences are of 
the order of the experimental precision. 

Yet, a closer look at the activation barriers and pref- 
actors (Tab. P), as derived from the experimental phase 
diagrams, poses a number of questions. First, why are 
the activation energies and prefactors in the adsorption 
and growth phase diagram so different, even though ab- 
solute values of transition fluxes are close? Second, why 
are the values for the same transition (2 — > 3 for growth) 
in Refs. ^| and ^3 so different? Finally, while the bar- 
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FIG. 7: Comparison between the experimental growth phase 
diagram in Ref. Qj] (open symbols) and the diagram derived 
from the the Ga adsorption data in this work (solid lines) . The 
different regimes 1, 2, and 3 of the adsorption phase diagram 
are indicated, which correspond to regimes B, C, and D in 
Ref. 0, respectively. 

rier of 2.8 eV for the 2^3 transition in the growth 
phase diagram in Ref. ^| is close to the cohesive energy 
of bulk Ga and has thus been directly interpreted as a 
Ga desorption barrier, what is the physical origin of the 
5.1 eV activation energy and the meaning of a prefactor 
of 10 25 Hz? The latter value is fundamentally different 
from prefactors typically observed/calculated for diffu- 
sion or desorption processes (which are of the order of 
- 10 13 Hz). 

A. Growth and adsorption model 

To address the above questions, we have analyzed the 
data in terms of a simple growth model and in combi- 
nation with first principles total energy calculations. In 
order to simplify the discussion of the growth model, we 
divide the problem in two parts: First, we derive how the 
density of Ga adatoms p on the surface^ depends on pa- 
rameters such as temperature T, Ga flux $, and growth 
rate v. Second, we calculate the critical adatom density, 
at which nucleation occurs and the system undergoes a 
phase transition. 

The adatom coverage is given by: 



adatoms and r gives the lateral position on the sur- 
face. For step flow growth (as realized in Ga-rich growth 
conditions!!), incorporation occurs essentially at the step 
edges. Since these are moving, the incorporation rate 
T inc ( r ' ^) * s m g enera l inhomogeneous and time depen- 
dent. For ideal step flow it will be zero on the terraces 
and > only at the step edges. The incorporation rate 
and the growth rate v are directly related by: 

"< f) = 3 /dro p{ '- t)dc ■ (5) 

Here the integration is performed over the total surface 
area A. Since in our experimental setup stationary con- 
ditions are realized, the explicit time dependence in the 
growth rate disappears (i.e. v(t) — v). Furthermore, in 
stationary conditions Eq. also simplifies: dp/dt = 0, 
i.e. the left hand side becomes zero. A further simplifica- 
tion in Eq. © can be made by taking into account that 
the phase transitions we are interested in occur exclu- 
sively under Ga-rich conditions. Under these conditions 
the surface steps on the surface are Ga-terminated 2 ^ and 
the incorporation rate Tj~* at such a step will be small. 
This is in contrast to the conventional step flow picture 
where the sticking probability of an atom at the step 
edge is assumed to be close to one. The difference to 
the conventional model is due to the fact that we have 
two species with very different concentrations. For ni- 
trogen (which is the minority species for very Ga-rich 
conditions) the sticking coefficient at the Ga-terminated 
step edges will be close to one and steps act as sinks 
to the nitrogen concentration which becomes highly in- 
homegenous. For Ga atoms, however, the sticking proba- 
blity is low. Thus, the effect of steps on Ga will be small 
and the Ga-adatom density is virtually homogeneous, i.e. 
DV 2 p(r) = 0. 

Based on the above discussion, Eq. ((3J) can be written 
as: 



= $ - v po ■ (6) 

Tides 

Here, p denotes the equilibrium adatom density. We 
note that this equation holds for both growth (v > 0) 
and adsorption (y = 0). The solution is easily found to 
be 



dp(r,t) 2 , , 1 , , 1 , . Po = (* - v) r dcs . (7) 

-W>,() + $ p(r,t) p{v,t) . 



dt 'Tine ^~dc 

(3) 



Nucleation occurs if the (stable) nuclei are in thermo- 
dynamic equilibrium with the lattice gas (which is de- 



Here, the desorption time T dcs is given by , u , , „ , ., r « , , 

' ^ scribed by the adatom density) on the surface. At low 

densities (po -C 1), interactions in the lattice gas itself 

_i _ -E dC s/k B T ( A \ can be neglected and we obtain 

T dcs — ^dese ° & , (4) 



with -Edcs the desorption barrier and kg the Boltzmann 
constant. D is the surface diffusion constant for Ga 



(8) 
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Here, N a< x is the number of adatoms in the lattice gas, 
N s ites 1S th e total number of surface sites which can be 
occupied by the adatoms, and AE nuc is the energy the 
adatom gains if it is attached to a subcritical nucleus 
making the latter stable. 

Based on the above model, we can directly obtain the 
critical Ga flux <& cr it at which the phase transitions occur. 
Combining Eqs. Q and JSJ) gives: 

This equation can be rewritten using Eq. (0J as 

<&c ri t=« + ^d Cse - (Bd - +A£ ™ )/fcBT (10) 

and applies both to the adsorption and growth phase 
diagram. The activation energy is thus expected to rep- 
resent the total binding energy of a Ga atom in a critical 
Ga cluster. 



B. First principles analysis 

A comparison with the experimental results (see 
Tab.QJ shows that energy and prefactor are not constant 
[as expected from Eq. IjlUI) ] but vary largely with the 
growth conditions. As a general trend, one finds that 
the activation energy and the prefactor decrease with in- 
creasing growth rate (N flux) . It is also interesting to note 
that only in the case of high growth rate (v = 1.1 ML/s) 
the prefactor (1 x 10 14 Hz) is close to the typical at- 
tempt frequencies observed/expected for desorption, i.e. 
in the 10 13 Hz range. For conditions where growth is 
slow (v ~ 0.28 ML/s) or absent (adsorption), prefactors 
are found which are many orders of magnitude larger. 

In order to identify the origin of these apparent dis- 
crepancies, we have explicitly calculated the desorption 
and the formation of small Ga clusters on the Ga bilayer 
surface employing density functional theory. In the fol- 
lowing we will focus on the Ga-bilayer structure (which 
corresponds to the 2 — > 3 transition). Based on the al- 
most identical energies of the 1^2 and 2-^3 tran- 
sitions in the adsorption phase diagram we expect the 
mechanisms/energetics to be rather similar. 

Specifically, we use soft Troullier-Martins 
pseudopotentials 24 and the Perdew-Burke-Ernzerhof 
generalized-gradient approximation (PBE-GGA) to de- 
scribe exchange/correlation^ The Ga 3c? semicore states 
were described in the frozen core approximation (nlcc)^& 
Details about the method can be found elsewhereiSI 
The calculations have been performed with a plane wave 
basis set (energy cutoff: 50 Ry). The Brillouin zone has 
been sampled by a 2 x 2 x 1 Monkhorst-Pack mesh for the 
(2V3 x 2\/3) unit cell and 4 x 4 x 1 for the (V5 x \/3) 
unit celli2& The Ga bilayer with one up to four adatoms 
has been modeled by a slab consisting of 2 double layers 
and (-\/3 x (to describe the adsorption of a single 




o#o 

# Ga adatom (layer 0) 

• Contracted Ga overlayer (layer 1 ) 
O Ga (layer 2) 

FIG. 8: Schematic top view of an adatom on the Ga bilayer 
structure. The white balls mark the positions of the Ga atoms 
in the second layer. The gray balls mark the positions of the 
Ga atoms in the contracted Ga epilayer and the black ball 
the Ga adatom in the T4 position. The dashed line shows the 
lxl surface unit cell of the ideal bulk truncated GaN (0001) 
surface. The solid line shows the surface unit cell of the Ga 
bilayer structure. 



adatom, see Fig. [SJ and larger (2-^/3 x 2\/3) (to describe 
adatom islands, see Fig. unit cells. Increasing the slab 
thickness to 4 double layers changes the surface energy 
by less than 0.5 meV per surface unit cell. Adatoms 
are added only on the upper surface of the slab. The 
lower surface has been passivated by pseudo-hydrogen 
to remove the electrically active surface states. The 
adatom(s) and the first two surface layers have been 
fully relaxed. Detailed convergence checks can be found 
in Refs. I2l3(ill3ll 




# Ga adatom (layer 0) 

• Contracted Ga overlayer (layer 1 ) 
O Ga (layer 2) 

FIG. 9: Schematic top view of a 2 atom island (dimer) on the 
Ga bilayer structure. The solid line shows the (2\/3 x 2\/3) 
surface unit cell which has been used to describe this struc- 
ture. 
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Based on these studies, we have calculated the des- 
orption energy of an adatom and the binding energy of 
adatoms in small islands. The desorption energy (the en- 
ergy needed to remove the adatom from the surface) is 
defined by 



/ Trradatom z^slab 77'atom\ 
^des - - (E tot - ^tot _ ^tot J 



(11) 



where E^ tonl is the total energy of the surface including 
the adatom, E^ \ h that of the free surface and E^° m that 
of the (spin-polarized 3 ^) Ga atom. Using this expres- 
sion we find an adatom binding energy of 2.52 eV for the 
(a/3 x \/3) structure and 2.41 eV for the (2a/3 x 2a/3) 
unit cell. In the equilibrium configuration, the adatom 
sits on a three-fold coordinated hollow site. The Ga-Ga 
bond length between adatom and surface layer is 2.68 A, 
i.e. very close to the nearest neighbor distance in a-Ga 
of 2.71 A^i 

The binding energy of an adatom in an island consist- 
ing of n a d adatoms is given by: 



isl 



1 

«ad 



/ ZTiisl TTTslab 77iatom\ 771 



(12) 



For an island consisting of 2 adatoms we find Ei S \ = 
0.15 eV. For larger islands consisting of 3 adatoms -Ej s i = 
0.30 eV and E isl = 0.34 eV for a 4 atom island. The 
numbers are the island formation energies as calculated 
in a (2\/3 x 2\/3) cell. It is interesting to note here 
that all islands are unstable against the formation of Ga- 
droplets: the formation energy of a Ga atom in an island 
{Edcs + Ei S \ ~ 2.75 eV) is smaller than the cohesive en- 
ergy of bulk Ga of 2.8 eVM Therefore the islands act as 
nucleation centers for Ga droplet formation. 



70 r 
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20 [ 



2.0 



3.0 



4.0 

-exp 



5.0 



6.0 



E7 (eV) 



FIG. 10: Natural logarithm of the experimental prefactors 
In (VfaH) as a function of the activation energy E^ p for various 
data in this work and Refs . Il4l and UtI (full circles; see Tab.^J. 
The solid line represents a linear fit. 



we assume a linear dependence: Ea = Eq + a(T — To) 
with Eq the temperature independent contribution, To 
the temperature offset, and a the linear temperature co- 
efficient. Equation l|10f> then becomes 



-a/k B p -{E -aT a )/k B T 



(14) 



A comparison between Eqs. (|13fl and (|14fl gives the fol- 
lowing relations: 



E^ p = Eq—o/Tq and 



This leads to 



M^def) = ln <Tdcs> 



a 

k B 
(15) 



C. Interpetation of the results 

We can now compare these energies with those ob- 
tained from the analysis of the experimental data 
(Tab.P). As can be seen, the activation energy is close to 
our calculated desorption energy only for the high growth 
case (l.lML/s). For lower growth rates or adsorption, 
activation energies are found which are way too large. A 
closer analysis of the prefactors and the activation ener- 
gies shows a clear relation (see Fig. I10|) : 

ln(C s P )=«i^ XP + «o • (13) 

The superscript "exp" indicates experimental data. For 
the values given in Tab. [I] we get ai = 11.45 eV -1 and 
a = -0.427. 

The observed relation between prefactor and activa- 
tion barrier can be explained in terms of a temperature 
dependent activation energy. Since the experimentally 
accessible temperature range is rather small (~ 50 K), 



a = m(^ dos ) - — — and d\ = — — . (16) 

Using these relations and Eq = 2.75 eV (as found from 
our first principles calculations), we obtain To = 740 °C 
(which is close to the experimentally accessed temper- 
ature range, so the linear approximation for the tem- 
perature dependence of Ea is well justified) and Vd C s = 
3x 10 13 Hz (which is close to typical attempt frequencies). 

Using the above parameters and Eq. I|15l) . the linear 
temperature coefficient a — (Eq — E a * p )/Tq can be com- 
puted for all transitions. The result is shown in Tab. [I] 
Based on these values we can renormalize the prefactors 
following z/g" = Vdcse~ a ^ kB ■ The resulting numbers cal- 
culated with ^d cs = 3x 10 13 Hz are listed in Tab. [I] The 
good agreement with the experimental prefactors i/LF 
shows that the model consistently describes all previous 
experimental studies. The large variation in frequencies 
and activation energies can be explained by assuming 
that the N flux affects the temperature dependence: it is 
largest if no N flux is present and monotonously decreases 
with increasing N flux (growth rate). 
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It is interesting to note that, although the temperature 
dependence has a huge effect on the experimentally mea- 
sured apparent activation barrier E c ^ v = E — aT (by 
almost a factor of two) and prefactor vf*? (by up to 12 
orders of magnitude), the actual change in the effective 
desorption energy Ea — Eq + a(T — To) is small: in the 
case of adsorption (v = 0), the activation barrier changes 
from 2.84 eV to 2.73 eV within the experimentally mea- 
sured temperature window (700 °C to 750° C). In the case 
of GaN growth, the variation is even smaller. This can 
be intuitively understood by considering that the exper- 
imentally measured activation energy represents a linear 
projection to zero temperature even when its real temper- 
ature dependence deviates strongly from linear behavior 
outside our experimental temperature window. 

V. CONCLUSION 

Based on a combination of experimental RHEED stud- 
ies and first-principle growth models we have (i) quanti- 
tatively determined the Ga coverage on the GaN (0001) 
surface during adsorption as a function of Ga flux and 
substrate temperature and (ii) derived a model, which 
consistently describes the adsorption of Ga on GaN sur- 
faces as well as the accumulation of Ga during Ga-rich 
GaN growth. This model resolves the discrepancy in pre- 
vious measurements of the activation energy characteriz- 
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